RBL2 represses the transcriptional activity of Multicilin to inhibit multiciliogenesis

A core pathophysiologic feature underlying many respiratory diseases is multiciliated cell dysfunction, leading to inadequate mucociliary clearance. Due to the prevalence and highly variable etiology of mucociliary dysfunction in respiratory diseases, it is critical to understand the mechanisms controlling multiciliogenesis that may be targeted to restore functional mucociliary clearance. Multicilin, in a complex with E2F4, is necessary and sufficient to drive multiciliogenesis in airway epithelia, however this does not apply to all cell types, nor does it occur evenly across all cells in the same cell population. In this study we further investigated how co-factors regulate the ability of Multicilin to drive multiciliogenesis. Combining data in mouse embryonic fibroblasts and human bronchial epithelial cells, we identify RBL2 as a repressor of the transcriptional activity of Multicilin. Knockdown of RBL2 in submerged cultures or phosphorylation of RBL2 in response to apical air exposure, in the presence of Multicilin, allows multiciliogenesis to progress. These data demonstrate a dynamic interaction between RBL2 and Multicilin that regulates the capacity of cells to differentiate and multiciliate. Identification of this mechanism has important implications for facilitating MCC differentiation in diseases with impaired mucociliary clearance.


INTRODUCTION
Multiciliated cells (MCCs) are key drivers of directional ciliary flow along epithelia in many tissues including the airways, reproductive tract, and brain [1][2][3][4][5][6].In the airways, MCCs play an integral role in the lung's primary innate defense system, providing both a physical barrier from the environment and mechanical force that drives mucus, pathogens, and foreign debris out of the lungs in a process called mucociliary clearance [7].Due to their central role in maintaining airway homeostasis, dysregulated MCC function is a critical pathophysiologic feature underlying many respiratory diseases, including chronic obstructive pulmonary disease (COPD) and cystic fibrosis (CF) [8,9].Primary ciliary dyskinesia (PCD) [10], a genetic disease arising from mutations in cilia related genes, results in non-motile, dyskinetic, or absent cilia, and symptomatically recapitulates secondary ciliary dyskinesia, characterized as the disruption of ciliary function from disease or injury [11][12][13][14].Repeated injury, characteristic of chronic lung disease, can lead to airway remodeling and changes in the cellular composition of airway epithelium, resulting in a decreased number of MCCs covering the airway surface [8].Due the prevalence of MCC dysfunction in respiratory diseases, it is therefore critical to understand the mechanisms controlling multiciliogenesis and how these may be targeted to restore functional mucociliary clearance.
To accomplish their function, each airway MCC must generate ~150 motile cilia [15].Cilia are microtubule-based organelles anchored to the plasma membrane by a basal body, formed from a mature parent centriole during cellular quiescence [15].Most cell types extend a single cilium because centriole number is tightly controlled, with centriole duplication occurring only once per cell cycle [16].By contrast, differentiating MCCs generate hundreds of centrioles as post-mitotic cells, thus engaging novel centriole biogenesis mechanisms outside the cell cycle and in a manner capable of significantly increasing centriole number.This process of massive centriole amplification is mediated by a MCC-specific transcriptional program that induces expression of the centriolar proteins that allow for centriolar duplication, which normally occurs only during S phase, as well proteins that form a structure unique to MCCs known as the deuterosome [17][18][19][20][21][22][23].
Geminin family proteins, GEMC1 and Multicilin, encoded by the gene MCIDAS, have been identified as master transcriptional regulators of multiciliogenesis, controlling both the induction of deuterosome-mediated centriole biogenesis and activating the gene regulatory network necessary to produce motile cilia [4,[24][25][26][27].GEMC1 and Multicilin share similar sequence homology, with a conserved coiled-coil domain and carboxy terminus; however, GEMC1 has been shown to help specify MCC fate, while Multicilin is sufficient to induce multiciliogenesis and can bypass GEMC1 [28].Upstream of GEMC1 and Multicilin, inhibition of NOTCH signaling is critical in regulating multiciliogenesis; a direct link between NOTCH-mediated cell fate determination and induction of Multicilin remains unknown [28].During multiciliogenesis, Multicilin acts transcriptionally in a complex with cell cycle transcription factors E2F4 or E2F5, along with dimerization partner DP-1, activating the expression of centriole biogenesis genes and key downstream transcription factors, including c-MYB, TP73, and FOXJ1 [24].Paradoxically, this activation occurs despite the well-known role of the E2F4 and 5 in post-mitotic cells as repressors of gene expression during cell cycle progression, by acting in a complex with the retinoblastoma (RB) family of corepressors, and the multi-vulval class B (DREAM) complex [29].
While Multicilin is known to be sufficient to drive multiciliogenesis in certain cell types [30][31][32], it has become apparent that this does not apply to all cell types, nor occur evenly across all cells in the same cell population.In our prior study using mouse embryonic fibroblasts (MEFs), expression of Multicilin alone was insufficient and could only stimulate multiciliogenesis after serum starvation and when co-expressed with E2f4 modified by replacing the C-terminus with a VP16 transactivator, indicating that the Multicilin requires a specific cellular context in order to promote multicilogenesis [33].In this report we further investigate how co-factors regulate the ability of Multicilin to drive multiciliogenesis.We identify RBL2 as a repressor of the transcriptional activity of Multicilin and demonstrate its dynamic interaction during differentiation of human airway epithelial cells at the airliquid interface (ALI).Identification of this mechanism has important implications for facilitating MCC differentiation in diseases where mucociliary clearance is impeded.

E2F4 C-terminus substitution with VP16 attenuates Rbl2-Multicilin interaction in MEFs
In their canonical roles, the E2F4 pocket protein binding domain (PPBD) is known to bind co-repressors to inhibit cell cycle gene expression, while the transactivating domain (TAD) of E2F family proteins associates with epigenetic modifiers to activate transcription (Supplementary Fig. S1a) [34].We have previously reported that overexpression of Multicilin was not sufficient to induce multiciliogenesis effectively in MEFs [33].This block could be overcome by overexpressing Multicilin along with a modified form of E2F4 where both the PPBD and TAD domains were replaced with a viral transactivation domain VP16 (E2f4ΔCT-VP16), as evidenced by robust centriole amplification and increased expression of multiciliogenesis-related proteins (Supplementary Fig. S1b, c) [33].To determine whether the proteins associated with the Multicilin complex change when cells express Multicilin along with wildtype E2F4 (limited multiciliogenesis) compared to E2f4ΔCT-VP16 (extensive multiciliogenesis) we performed proteomics.
Adenoviral vectors were used to express Multicilin tagged at the N-terminus with 3xFLAG (FLAG-MCI, Fig. 1a), FLAG-MCI along with wildtype E2F4 (FLAG-MCI+E2f4, Fig. 1b), or FLAG-MCI along with E2f4ΔCT-VP16 (FLAG-MCI+E2f4VP16, Fig. 1c).MEFs were transduced, and immunoprecipitation performed using FLAG to copurify proteins complexed to Multicilin.The bound proteins to FLAG-MCI were eluted and silver staining performed to identify the MCI-3xFLAG (Fig. 1e).Mass spectrometry analysis identified multiple proteins with similar abundance between FLAG-MCI and FLAG-MCI+E2f4VP16, with 8 of the 10 most abundant proteins being present in both conditions.These proteins included known Multicilin binding partners, TFDP-1, GMNN, E2F4, and E2F5, as well as ribosomal proteins whose interaction with Multicilin has not previously been described: RPS27, RPS27L and RPS14 (Fig. 1d).A complete list of the proteins complexed with Multicilin is included in Supplementary Table S1.The other 2 of the 10 most abundant proteins bound to FLAG-MCI were 26 S protease regulatory subunit 7 (PSMC2) and Retinoblastoma-like protein 2 (RBL2).Notably, RBL2 had the largest decrease in relative abundance comparing FLAG-MCI to E2f4VP16 of all proteins detected in both samples, with a ratio of 9.2:1.
Immunoblotting analysis detected a strong interaction between Multicilin and RBL2 in MEFs expressing either FLAG-MCI or FLAG-MCI+E2f4; however, the RBL2-Multicilin interaction was abolished in MEFs expressing FLAG-MCI+E2f4VP16 (Fig. 1f).Similarly, an  S1). e Representative silver stain of total protein lysate, the red arrowhead indicates 3xFLAG-MCI.f Representative western blot analysis for Rbl1, Rbl2, and FLAG for total protein (top) and co-immunoprecipitations (Co-IP, bottom) with FLAG from non-infected mouse embryonic fibroblasts (MEFs) and MEFs infected with Ad5 FLAG-MCI, FLAG-MCI+E2f4 and FLAG-MCI+E2f4VP16.
interaction with RBL1 was observed in MEFs expressing FLAG-MCI in the presence of E2f4, but not E2f4VP16 (Fig. 1f).The retinoblastoma family of proteins, including RB, RBL1 and RBL2, have been well characterized in their roles as cell cycle inhibitors through interactions with the E2F family of transcription factors [35,36].The RBL2-E2F4 complexes have been identified as the predominant repressors of cell cycle genes during quiescence [35,36].Taken together, the canonical role of RBL2 as a corepressor of E2F4 and our observation that RBL2 complexes strongly with Multicilin in the presence of E2f4 relative to E2f4ΔCT-VP16, suggest that RBL2 may regulate the capacity of the Multicilin transcriptional complex to activate MCC gene expression.

RBL2 represses the ability of Multicilin to drive ectopic multiciliogenesis in MEFs
To further investigate a regulatory role of RBL2 in the Multicilin-E2F4 transcriptional complex, we tested siRNAs for Rbl2 knockdown in MEFs (siRbl2#1-3, Supplementary Table S2), identifying siRbl2#1 (siRbl2) as most efficient at reducing Rbl2 gene expression, with a reduction of 1.42 ± 0.44 log2 fold change compared to a nontargeting control (siCTL) (Supplementary Fig. S2a).An siRNA targeting Rbl1 (siRbl1) was also tested and compared as an additional control.Successful knockdown was validated at the protein level by western blot (Fig. 2a, Supplementary Fig. S2b, c).
To further evaluate the ability of RBL2 to regulate Multicilindependent activation of multiciliogenesis, we performed bulk RNA-sequencing to determine Multicilin-induced transcriptional changes comparing Multicilin induced MEFS with siCTL, siRbl2 or E2f4VP16 to noninfected MEFs (Supplementary Table S3).Multicilin-induction with E2f4VP16 induced the most dramatic changes in overall gene expression compared to the non-infected MEFs, with similar trends occurring in the presence of siRbl2 and with siCTL having less substantial transcriptional changes (Supplementary Fig. S3).Evaluating a panel of MCC-related genes [37], both E2f4VP16 and siRbl2 induce comparable transcriptional changes, while siCTL more closely reflects non-infected MEFs (Fig. 3a, Supplementary Fig. S4).Gene Ontogeny (GO) analysis of differentially expressed genes (DEGs) comparing experimental conditions and the non-infected controls (adj p > 0.05, |FC | > 2) highlighted significant enrichment in Multicilin-induced MEFs, in the presence of either E2f4VP16 or siRbl2, for biological processes (GO:BPs) such as cilium organization, cilium movement, and protein localization to cilium as well cellular components (GO:CCs) cilium, centriole, manchette and dynein axonemal particle (Fig. 3b).Multicilin-induction with siCTL had a significant enrichment in GO:CC deuterosome and GO:BPs cilium organization and multiciliated differentiation; however, the most significantly enriched GO:BPs were immune system process and response to interferonbeta (Fig. 3b).Similar numbers of DEGs from these GO terms were also noted in the E2f4VP16 and siRbl2 MEFs, however the greater total number of DEGs decreased the enrichment significance for these GO terms (Supplementary Table S3).Comparing common DEGs between samples, Multicilin-induced MEFs with E2f4VP16 have higher MCC gene expression compared to Rbl2-depleted MEFs (Fig. 3c), however both are significantly increased compared to siCTL (Fig. 3d, e).The Venn diagram shown in Fig. 3f compares all DEGs, with >2-fold changes in expression, across all experimental conditions relative to non-infected controls.We identified 448 DEGs significantly upregulated specifically in the Multicilininduced and Rbl2-depleted or E2f4VP16 expressing MEFs, but not in the Multicilin-induced siCTL MEFs (Fig. 3g).Of these 448 genes, 99 are present in a previously reported mouse MCC-associated gene list (Fig. 3g, Supplementary Table S4) [37,38].This data supports a regulatory role for the suppression of Rbl2 in enabling Multicilin to activate transcriptional pathways to induce massive centriole biogenesis and form multiple motile cilia.

RBL2 expression does not change the extent of multiciliogenesis from HBECs
Since knockdown of Rbl2 allowed Multicilin to induce ectopic multiciliogenesis in MEFs, a cell type that is not innately capable of multiciliogenesis, we investigated whether RBL2 knockdown in human MCC progenitors can potentiate multiciliogenesis.To accomplish this, we compared the differentiation of human bronchial epithelial cells (HBECs) at the air-liquid interface (ALI) with or without RBL2 knockdown (Fig. 4).qRT-PCR analysis after 14 days of ALI differentiation confirmed that shRBL2 (shRBL2#1 and #2) significantly decreased expression of RBL2 when compared to non-targeting controls (shNTC#1 and #2), by 2.37 ± 0.63 and 1.30 ± 0.19 log2 fold change, respectively (Fig. 4a).The reduction in RBL2 expression resulted in a modest increase in the expression of early MCC differentiation genes including SAS6, DEUP1 and TP73 (Fig. 4b), however the expression of genes associated with differentiated human airway epithelium, including goblet (MUC5AC), club (CC10), and ciliated (FOXJ1) cells was unchanged (Fig. 4c).This suggests that the proportion of MCCs generated in the presence of RBL2 knockdown was unchanged.Immunoblotting confirmed knockdown of RBL2 at the protein level alongside increased expression of TP73, while more mature MCC markers, FOXJ1 and axonemal protein CCDC39, remained unchanged (Fig. 4d, e).Immunostaining for cilia maker acetylated α-tubulin (ATUB) was used to evaluate the efficiency of MCC differentiation by quantifying the area of ciliated cell coverage on the apical surface at ALI day 14, while live video analysis was used to measure cilia beat frequency (CBF).No significant changes in total ciliated cell coverage, nor in CBF were detected in RBL2 knockdowns compared to controls (Fig. 4g, h).These data show that ectopic reduction of RBL2 expression does not change MCC maturation and function and suggests that RBL2 activity may already be endogenously regulated during the early stages of HBEC differentiation at the ALI.

ALI culture phosphorylates RBL2, allowing Multicilin to drive multiciliogenesis in HBECs
To determine whether RBL2 interactions with Multicilin are indeed altered during the early phase of ALI differentiation of HBECs, we evaluated RBL2-Multicilin interactions comparing ALI to submerged cultures in differentiation media.Doxycycline (DOX) induction of lentiviral-transduced FLAG-tagged Multicilin in HBECs was performed over a 4-day period prior to sample collection (Fig. 5a, Supplementary Fig. S5a-b).Immunoprecipitation of FLAGtagged Multicilin and immunoblotting for RBL2 revealed a higher amount of RBL2 complexed with Multicilin in the submerged cultures (Fig. 5a, Sub D0 and Sub D4), compared to ALI cultures (Fig. 5a, ALI D4 and ALI D21).Total RBL2 remained consistent in all conditions (Fig. 5a, Supplementary Fig. S5b).These data demonstrate that ALI culture attenuates the direct interaction between RBL2 and Multicilin in HBECs, which in turn potentiates Multicilin transcriptional activity to drive MCC differentiation.
Cyclin-dependent kinases are known to phosphorylate the RB family of proteins, interrupting their interaction with E2F family transcription factors to regulate gene expression and cell cycle progression [39,40].Canonically, RBL1/2 interactions with E2F4/ 5 actively repress cell cycle genes [36].The ability of Multicilin to co-opt E2F4 to induce centriole biogenesis in cells not actively dividing, suggests a mechanism whereby RBL2 may be phosphorylated and inactivated.RBL2 has over 26 reported phosphorylation sites [41][42][43], of which hyperphosphorylation is associated with dissociation from E2F4 and RBL2 degradation [41,44,45].Among these phosphorylation sites, phosphorylation of S672 has been shown to be important for attenuating E2F4-RBL2 interactions and is necessary for ubiquitination and degradation [44].We therefore investigated whether S672 RBL2 phosphorylation status differed between HBECs in ALI and submerged culture.We observed no significant difference in total RBL2 protein expression while S672-phosphorylated RBL2 (S672-RBL2) was significantly higher in HBECs cultures at the ALI compared to submerged (Fig. 5b-d).Immunoblotting analysis   S3) comparing non-infected MEFs and FLAG-MCI MEFs with either siCTL, siRbl2 or E2f4VP16.GO terms include biological process (GO:BP) and cellular components (GO:CC).c-e Scatter plots comparing significant DEGs in FLAG-MCI MEFs (relative to uninfected, absolute log2 fold change>2, adj.p < 0.5) for siRbl2 vs. siCTL (c), E2f4VP16 vs siCTL (d), and siRbl2 vs E2f4VP16 (e).f, g Venn diagrams representing significantly increased DEGs in FLAG-MCI relative to uninfected controls (log2 fold change>2, adj.p < 0.5, Supplementary Table S5) comparing siRbl2, siCTL and E2f4VP16 (f) and compared to known MCC-specific genes [37] (g).N = 3 experimental replicates for all RNAseq data.confirmed that MCC differentiation was activated, evidenced by the induction of Multicilin transcriptional targets, FOXJ1 and TP73, at D4 ALI differentiation (Fig. 5b).These data suggest that the RBL2-dependent regulation of Multicilin activity, during ALI differentiation, is associated with increased RBL2 phosphorylation.
Expression of Multicilin in is known to be dependent on the inhibition of NOTCH signaling downstream of air exposure [13].To determine if RBL2 phosphorylation in differentiating HBECs can be regulated by NOTCH inhibition and Multicilin activity, independent of ALI, we inhibited NOTCH signaling using 5 µM of the γ-secretase inhibitor (DAPT) and induced Multicilin for 4 days in HBECs submerged in differentiation media.Immunoblotting analysis revealed that inhibition of NOTCH signaling alone decreased total RBL2 protein expression, while Multicilin induction alone increased S672-RBL2 expression; both conditions resulted in increases in S672-RBL2:total RBL2 (Fig. 5e, f).Interestingly, only the combination of NOTCH-inhibition and Multicilin-induction significantly increased the S672-RBL2:total RBL2 compared to each factor alone (Fig. 5f, g).This data supports a regulation of RBL2 phosphorylation downstream of both NOTCH-inhibition and Multicilin-induction and suggests that they work synergistically to phosphorylate RBL2 at S672.

RBL2 knockdown increases the transcriptional activity of Multicilin in submerged HBECs
As RBL2 phosphorylation is associated with dissociation from the E2F proteins allowing transcriptional activation, we depleted RBL2 in submerged Multicilin-inducible HBECs to determine whether Multicilin-dependent transcriptional activation of multiciliogenesis would be potentiated independent of regulation by ALI differentiation.Multicilin expression was induced for 4 days in submerged cultures.Induction alone (red bars, Fig. 6a) had no effect on RBL2 expression, but significantly increased expression of MCC genes TP73, endogenous MCIDAS, and FOXJ1 (Fig. 6a).RBL2 knockdown alone had no significant impact on TP73 gene or protein expression (Fig. 6a-d), however, increased MCC gene expression was observed for FOXJ1 (Fig. 6a).Interestingly, the combination of RBL2-knockdown and Multicilin-induction led to a significant increase in the expression of MCC genes SAS6, DEUP1, TP73 and FOXJ1 (blue bars + DOX, Fig. 6a), the number of TP73positive cells (Fig. 6b, c), and total TP73 protein expression (Fig. 6d) when compared to either condition alone.Additionally, only in the combination of RBL2-knockdown and Multicilin-induction were we able to observe some cells with multi-nucleated basal bodies and visible cilia axonemes, shown in confocal images for pericentrin and ATUB (Fig. 6e).We also found that Multicilininduced MCCs with RBL2 depletion had active ciliary beating, observed in Supplementary Movie S1, indicating functional cilia differentiation.These data support a new role for RBL2 in the direct inhibition of the activity of Multicilin.Only in the absence of RBL2 can Multicilin efficiently drive the transcriptional network necessary for multiciliogenesis.

DISCUSSION
This study was designed to elucidate the mechanisms that regulate the potential of human airway epithelial cells to become multiciliated.We previously reported that Multicilin can drive ectopic multiciliogenesis in MEFs, but only when using a modified version of E2f4 where the C-terminus is replaced with a VP16 transcriptional activator [33].Here we show that E2f4-VP16 disrupts the interaction of RBL2 with the Multicilin/E2F protein complex, and that depletion of RBL2 alone increases the transcriptional activity of Multicilin in both MEFs and HBECs.Additionally, culture of HBECs at the ALI increases the levels of phosphorylated S672-RBL2 and disrupts the interaction of Multicilin with RBL2.This phenotype was recapitulated in submerged cells only in the presence of both Multicilin-induction and NOTCHinhibition.Collectively, the data support a novel mechanism whereby the ALI facilitates phosphorylation of RBL2, downstream of the inhibition of NOTCH signaling, allowing the Multicilin-E2f4 transcriptional complex to drive multiciliogenesis independently of cell-cycle mediated RBL2 phosphorylation.This new mechanism is summarized in the schematic in Fig. 7.
Multicilin and GEMC1 initially were associated with the cell cycle as Geminin superfamily proteins through sequence homology of their coiled-coiled domains with Geminin, a negative regulator of DNA replication [24,46].Data from a variety species, however, indicates that Multicilin and GEMC1 expression is necessary for proper MCC differentiation [17,[24][25][26][27][46][47][48][49], with loss of function in either gene causing reduced generation of motile cilia [27,49], without apparently leading to general cell cycle defects.Both GEMC1 and Multicilin lack a known DNA binding domain, and instead use their conserved C-terminal TIRT domains to interact with the transcription factors E2F4 and E2F5 [28,49].GEMC1 has been shown to more effectively complex with E2F5 and is thought to specify MCC fate by driving the expression of Multicilin and other downstream MCC transcriptional regulators, while Multicilin complexes with both E2F4 and E2F5 [28,49] to potentiate multiciliogenesis through a positive feedback loop, activation of massive centriole biogenesis, and upregulating the expression of downstream transcription factors.E2F4 and E2F5 have been well characterized as repressors of cell cycle progression, through their interactions with RBL1 and RBL2, which recruit various epigenetic regulators to repress the expression of cell cycle genes [50][51][52][53].The mechanism by which Multicilin co-opts E2F4 to induce transcription is largely unexplored, however recent work has shown that the GEMC1 and MCIDAS complex with TP73 [54] and TRRAP [49] and that their TIRT domains confer differential affinity for SWI/SNF subcomplexes [49].E2F4 and E2F5 are dualrole transcription factors, both acting in essentially all cell types as repressors during the cell cycle, but as activators in progenitors during multiciliogenesis; how the cell differentially regulates the function of E2F4 and E2F5 in this context is still unknown.
Our previously published work demonstrated that multiciliogenesis can be induced in MEFs, a mesenchymal cell type that is not innately capable of MCC differentiation, through a process requiring serum starvation and expression of Multicilin and E2f4ΔCT-VP16 [33].This suggested that a specific cellular state is required to potentiate multiciliogenesis downstream of Multicilin induction, specifically requiring G0/G1 arrest as well as additional unknown factors that can facilitate the role of E2f4 as a transcriptional activator.Here our identification of Rbl2 as a repressor of Multicilin, an interaction which can be interrupted by replacing the PBDB with VP16, is not surprising as serum starvation has been shown to increase hypo-phosphorylation of Rbl2 and E2f-Rbl2 interactions in the MEF derived cell line, NIH3T3 [55].As Rbl2-depleted MEFs expressing FLAG-MCI were able to generate cilia akin to that of Multicilin-induced MEFs expressing E2f4VP16, it is likely that MEFs have an innate capacity to utilize endogenous E2f4 to function as a transcriptional co-activator, however, this is reduced by the presence of active Rbl2.
Interestingly, Rbl1 was not efficiently co-immunoprecipitated with FLAG-MCI when expressed alone in MEFs.Multicilin-RBL1 interactions did not appear to be altered by ALI differentiation (data not shown), despite RBL1, and not RBL2, being identified in complex with Multicilin and GEMC1 in recently reported experiments performed in HEK293 cells [49].RBL2 has been reported to associate with E2F4 during G0 and G1, while RBL1/E2F4 complexes form mainly during S-phase of the cell cycle [56].This may explain why we detected Rbl1 complexed with Multicilin only when wildtype E2f4 was overexpressed, as endogenous E2f4 would preferentially bind Rbl2, but in the presence of excess E2f4, complexes with E2f4 were able to form.Additionally, our Rbl1 knockdown did not increase MCC specification.Taken together this points to a mechanism where RBL1 regulates GEMC1/ E2F5 specification of MCCs, while RBL2 regulates the potential of Multicilin/E2F4 to drive multiciliogenesis.
In conclusion, we provide evidence to support a mechanism by which ALI culture decreases the interaction of RBL2 with Multicilin.Submersion is believed to mimic a hypoxic environment and inhibit differentiation by maintaining NOTCH signaling [13].We show here that RBL2 interaction with Multicilin is regulated by its phosphorylation and only the combination of NOTCH-inhibition and Multicilin-induction was capable of increasing phospho-RBL2 in submerged cultures, suggesting that RBL2 phosphorylation is downstream of NOTCH-inhibition.

MATERIALS AND METHODS Cell culture
Mouse embryonic fibroblasts (MEFs).Freshly prepared primary MEFs were obtained from the Genome Manipulation Core at the Salk Institute (strain DR4, non-irradiated).MEFs were grown in DMEM supplemented with 10% FBS and 0.1 mM nonessential amino acid (#11140-050; Invitrogen) in a 37 °C humidified incubator with 5% CO2.Only MEFs cultured for less than 4 passages were used in experiments.
Human bronchial epithelial cells (HBECs).Primary HBECs were either a gift from Dr. Steve Brody (Wash U) or isolated from explant donor lungs as previously described [57] and cultured for one passage in Airway Epithelial Growth Medium (Promocell C-21160).Continued passaging of HBECs was performed in Pneumacult-Ex Plus (StemCell Technologies #05040) and differentiation of cells was performed in differentiation media, either PnuemaCult-ALI (StemCell Technologies #05001) or VALI media [58] on Transwell inserts (Corning #3470, #3460).Expansion and differentiation of HBECs was performed according to the Pneumacult-ALI manufacturer guidelines, with the exception that all cell culture surfaces were coated in PureCol (Advanced Biomatrix #5005) according to product guidelines.The use of de-identified HBECs was approved by the IRB of both the University of Southern California and the University of Iowa as "not Human Subjects Research."

Virus production and infection
pAd/CMV/V5 expressing mouse Multicilin and E2f4 proteins.The plasmids used in this study are modified from those described in Kim et al. [33].Mouse Multicilin (GenBank: AK134107.1)was obtained from the Riken mouse FANTOM clone library (clone ID 5830438C23).Mouse E2f4 (NCBI Ref: NM_148952.1)was derived from a cDNA obtained from Dharmacon (clone ID 4987691), by using the full-length coding domain (wildtype) or by using a C-terminal truncation (1-774 bp) that was replaced with the transcriptional activation domain of VP16 (Viral Protein 16; amino acids 413-490) from the UL48 gene of Herpes Simplex Virus-1 (HSV-1; GenBank: KM222726.1).The NLS based on SV40 T-antigen (PKKKRKV) and the 6xmyc tags were derived from the CS2 vectors.S2 cleavage sequence used in this study was T2A.Genes assembled in a pENTR/D-TOPO vector (K240020, Invitrogen) were validated by sequencing, and then transferred to the pAd/CMV/V5-DEST vector (#V49320; Invitrogen) using Gateway cloning.
Adenovirus.Adenoviruses were generated by GT3 core at the Salk Institute.Vector DNAs were initially transfected into 293 T cells, validated for intact protein expression by immunoblot and immunostaining, and then used to generate crude adenovirus lysates and titer calculated.Primary MEFs were plated with 10% FBS-DMEM for overnight, then treated with 2% FBS-DMEM for 1 day before adenovirus infection.Adenovirus infection was performed by adding the adenovirus crude lysate to the MEF cells for four hours, washing with pre-warmed PBS once, then growing infected cells with 2% FBS-DMEM.The cells were subjected to relevant analysis according to days post-infection (PI).
RBL2 knockdown via shRNA transduction.HBECs cultured in Pneumacult-Ex Plus were dissociated into single-cell suspensions using Accutase (#AT 104, Innovative Cell Technologies), counted, and the appropriate number of cells per infection condition were pelleted by centrifugation at 400 g for 5 minutes.Cells were then resuspended at a density of 500 K cells/mL in lentivirus diluted to an MOI of 20 in Pneumacult-Ex Plus with 10 µg/mL polybrene (#TR-1003, Millipore Sigma) and incubated at 37 °C for 2 hours, vortexed every 20 minutes to prevent cell clumping.Without removing the virus, cells were seeded at the appropriate cell density for either cell expansion or air-liquid interface differentiation and then incubated overnight 37 °C.Virus was removed the next day and cells were cultured in 1 µg/mL puromycin for at least three days before use in experiments.

Real-time quantitative polymerase chain reaction (qRT-PCR)
Total RNA was isolated from experimental MEFs (Quick-RNA prep kit; ZYMO RESEARCH, #R1054) or HBECs (Quick-RNA Microprep Kit; ZYMO RESEARCH, #R1051) and at least 500 ng RNA was converted in cDNA using High-Capacity cDNA Reverse Transcription Kit (Applied Biosystems, # 4368814).Gene expression was then assayed by quantitative PCR using PowerUp SYBR (Applied Biosystems, #A25742) and appropriate primer pairs in triplicate on the ABI Prism 7900HT thermal cycler (Applied Biosystems, #4329001) or the QuantStudio5 (Applied Biosystems, # A28575).Gene expression in HBECs and MEFs was normalized using reference genes RPLP0 and Gapdh, respectively and calculated as log2 fold change versus the average of control samples as stated.Primer sequences are listed in Supplementary Table S2.

Immunoprecipitation, immunoblotting analysis, and silver staining
For identifying proteins interacting with FLAG-tagged Multicilin by MASS analysis, MEFs were infected by pAd/CMV/V5 vector encoding 3xFLAG-Multicilin, NLS-6xmyc-E2f4ΔCt-VP16-T2A-Multicilin, and NLS-6xmyc-E2f4WT-T2A-Multicilin.Non-infected MEFs were used as the negative control.Two days after infection, the noninfected MEFs or infected MEFs were lysed in lysis buffer (50 mM HEPES-KOH [pH 7.5], 500NaCl, 5 mM EDTA [pH 8.0]), 2% Triton, DNase I (10ug/ml)) by incubation on ice for 20 min.Lysates were cleared by centrifugation at 12,000 rpm for 20 min at 4 °C, then precleared by incubation with Protein G agarose (Invitrogen #20397) agitating for 1 h at 4 °C.Spun-down beads were incubated with FLAG M2 agarose (Sigma #A2220) agitating for 2.5 hours at 4 °C.The beads were washed twice with wash buffer (1:1 mix of 50 mM HEPES and Lysis buffer), washed with lysis buffer twice, then washed twice again with wash buffer.Spun-down beads were incubated with 40uL of FLAG peptide (500 ng/mL) for bound protein elution overnight at 4 °C.Elutes were subjected to protein sampling for further immunoblot analysis and silver staining.Eluted protein samples and experimental protein lysates were subjected to SDS-PAGE and transferred to a PVDF membrane that was then blocked with 0.5% casein blocker (#161-0783, Bio-Rad) in PBS for 30 min followed by incubation with the indicated primary antibodies overnight at 4 °C in 0.1% Tween-20 in blocking buffer.After extensive washing in 0.1% Tween-20 in PBS, the blots were incubated with Alexa 680 or 800-conjugated anti-mouse or rabbit secondary antibodies (Invitrogen, A21058, A21076, A32735) for 45 min at room temperature, washed with 0.1% Tween-20 in PBS, and imaged using Odyssey (LI-COR).Silver staining was performed with Pierce Silver Stain Kit (Thermo #24612) based on the manufacturer's instructions.

Mass spectrometry analysis
Samples were precipitated by methanol/ chloroform and redissolved in 8 M urea/100 mM TEAB, pH 8.5.Proteins were reduced with 5 mM tris(2carboxyethyl)phosphine hydrochloride (TCEP, Sigma-Aldrich) and alkylated with 10 mM chloroacetamide (Sigma-Aldrich).Proteins were digested overnight at 37 °C in 2 M urea/100 mM TEAB, pH 8.5, with trypsin (Promega).Digestion was quenched with formic acid, 5% final concentration.The digested samples were analyzed on a Fusion Orbitrap tribrid mass spectrometer (Thermo).The digest was injected directly onto a 30 cm, 75 um ID column packed with BEH 1.7um C18 resin (Waters).Samples were separated at a flow rate of 300 nl/min on a nLC 1000 (Thermo).Buffer A and B were 0.1% formic acid in water and 0.1% formic acid in 90% acetonitrile, respectively.A gradient of 1-35% B over 110 min, an increase to 50% B over 10 min, an increase to 90% B over 10 min and held at 90%B for a final 10 min was used for 140 min total run time.Column was reequilibrated with 15 µl of buffer A prior to the injection of sample.Peptides were eluted directly from the tip of the column and nanosprayed directly into the mass spectrometer by application of 2.5 kV voltage at the back of the column.The Orbitrap Fusion was operated in a data-dependent mode.Full MS scans were collected in the Orbitrap at 120 K resolution with a mass range of 400 to 1600 m/z and an AGC target of 5e5.The cycle time was set to 3 sec, and within these 3 s the most abundant ions per scan were selected for CID MS/MS in the ion trap with an AGC target of 1e4 and minimum intensity of 5000.Maximum fill times were set to 50 ms and 100 ms for MS and MS/MS scans, respectively.Quadrupole isolation at 1.6 m/z was used, monoisotopic precursor selection was enabled, and dynamic exclusion was used with an exclusion duration of 5 sec.Protein and peptide identification were done with Integrated Proteomics Pipeline -IP2 (Integrated Proteomics Applications).Tandem mass spectra were extracted from raw files using RawConverter [62] and searched with ProLuCID [63] against Uniprot mouse database.The search space included all fully tryptic and half-tryptic peptide candidates.Carbamidomethylation on cysteine was considered as a static modification.Data was searched with 50 ppm precursor ion tolerance and 600 ppm fragment ion tolerance.Identified proteins were filtered to using DTASelect [64] and utilizing a target-decoy database search strategy to control the false discovery rate to 1% at the protein level [65].E.J. Quiroz et al.

Fig. 5
Fig. 5 Air-liquid interface attenuates the interaction between RBL2-Multicilin.a Representative western blot for total protein and immunoprecipitation with FLAG in FLAG-MCI induced HBECs either submerged (Sub) or ALI culture.b Representative western blot for total RBL2, phosphorylated serine 672 RBL2 (RBL2 S672) and MCC proteins TP73 and FOXJ1 at D0, D4 (Sub) and ALI D4.Protein quantification for total-RBL2 (c) and RBL2 S672:total RBL2 (d), N = 6 biological replicates.e Representative western blot for Multicilin-induced HBECs at Sub D4 in the presence of 5 µM DAPT or vehicle.f Protein quantification for DAPT treatment from N = 3, from 1-2 biological replicates and 1-2 experimental replicates.Data represents mean ± SEM and is compared using a paired t-test (c, d) and 2-way ANOVA and post-hoc Tukey's multiple comparisons test (f).

Fig. 7
Fig. 7 Regulation of Multiciliogenesis and the activity of Multicilin by RBL2.Schematic depicting the new mechanism of regulating MCI activity during multiciliogenesis.Unphosphorylated RBL2 binds to E2F4 and inhibits Multicilin-induction of MCC gene expression.Differentiation at an air-liquid interface and inhibition of NOTCH signaling both phosphorylate RBL2 and dissociate it from the E2F4/Multicilin complex, allowing for induction of MCC gene expression.